The Influence of Interference on the Kondo Effect in a Quantum Dot 
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We study the Kondo effect in a model system of a quantum dot embedded in an Aharanov- 
Bohm ring connected to two leads. By transforming to the scattering basis of the direct inter- 
lead tunneling, we are able to describe precisely how the Kondo screening of the dot spin occurs. 
We calculate the Kondo temperature and zero-temperature conductance and find that both are 
influenced by the Aharanov-Bohm ring as well as the electron density in the leads. We also calculate 
the form of an additional potential scattering term that arises at low energies due to the breaking of 
particle-hole symmetry. Many of our results are supported by numerical analysis using the numerical 
renormalization group. 



I. INTRODUCTION 

A quantum dot in a gated semi-conductor heterostruc- 
ture separating ballistic leads is known to exhibit some 
remarkable phenomena. Most strikingly, at low temper- 
atures, the conductance through the quantum dot in- 
creases as the temperature is lowered and can reach the 
ideal value of 2e 2 /h when the dot contains an odd number 
of electrons^—. This is due to the Kondo effect, involv- 
ing the screening of the spin 1/2 of the quantum dot by 
conduction electrons in the leads. If there is an addi- 
tional tunneling path connecting the two leads that does 
not pass through the quantum dot, then some interesting 
interference phenomena take placed—. 

Previous theoretical work on this problem has stud- 
ied both extended Aharonov-Bohm rings as well as short 
"Kondo-Fano" devices. The conductance^^ and ther- 
mopower— was found to exhibit an asymmetric Fano- 
like dependence on the energy level of the quantum dot. 
When the quantum dot is tuned to the Kondo regime that 
favors a local moment, a flux dependent Kondo tempera- 
ture has been proposed using different methods^— and 
the high and low-temperature conductance has been de- 
scribe d 10 ' 12 ^ . While some Numerical Renormalization 
Group (NRG) work was reported, this only studied the 
electron occupancy of the quantum dot 10 or the density 
of states on the quantum dot—, both of which can be 
approximately related to the conductance. It should also 
be noted that most of these studies assume a particle- 
hole symmetric dispersion relation and Fermi energy in 
the leads. 

In this paper, we reexamine the Kondo-Fano device 
using a combination of analytic and NRG methods. We 
only consider the Kondo regime where a local moment is 
favored on the quantum dot. We are able to reproduce 
many of the published results cited above as well as pre- 
dicting for the first time non-trivial dependence of the 
Kondo temperature and zero-temperature conductance 
on the electron density in the leads. Such a dependence 
on electron density has not been investigated before given 
that a particle-hole symmetric Fermi energy has always 
been assumed in the leads. We calculate the generation 
of additional potential scattering terms that have often 
been neglected in previous studies but which do lead to 



small corrections to the zero-temperature conductance. 
Numerical confirmation of many of our results is pro- 
vided for the first time using the NRG. 

Our analytic approach begins in § [TTJ with a tight- 
binding version of the Anderson model together with a 
direct tunneling term between the two leads and factors 
representing magnetic flux between the two conducting 
paths. Following refs. [lj, Il9l - [2ll we then perform an ex- 
act transformation to the "scattering basis" which diago- 
nalizes the direct tunelling part of the Hamiltonian when 
the hybridization to the Anderson impurity is turned off. 
This gives a Hamiltonian containing no direct tunnel- 
ing term, only the hybridization to the impurity, albeit 
with a more complicated dependence on flux, inter-lead 
tunneling, and particle momentum. The initial Hamilto- 
nian contains two scattering channels, the even and odd 
states, for example. After transforming to the scattering 
basis, only one linear combination of these appears in the 
Anderson hybridization; we refer to it as the "screening 
channel" . 

As we are primarily interested in the Kondo regime of 
the quantum dot, we perform a Schrieffer- Wolff trans- 
formation in the screening channel basis to obtain an 
effective Kondo model with an additional potential scat- 
tering term Kr that is of order the bare Kondo coupling 
and which vanishes (to this order) when the dot level is 
tuned to the symmetric value of Ed = — U/2 (these terms 
are defined in eq. <j5j) ) . This latter term is discussed in 
§ IIII Al Both the generated Kondo interaction and the 
potential scattering depend on the flux ip, the strength 
of the direct inter- lead coupling t', and the momentum of 
electrons in the leads. From the strength of this Kondo 
interaction we are able to obtain the dependence of the 
Kondo temperature on these model parameters as dis- 
cussed in § IIV Al 

Next, in § 1111 Bl we integrate out high energy states to 
obtain a low energy effective Hamiltonian. In addition to 
renormalizing the Kondo interaction, this also generates 
a small potential scattering term, Vr, of second order 
in the bare Kondo coupling. Hence, Vr contributes to 
the leading order term in the potential scattering when 
Ed = — U/2 though there may be other contributions as 
we discuss in the text. Otherwise, it is the Kr term dis- 
cussed above that provides the leading order contribution 



to the potential scattering. 

In § II V Bl we calculate the low temperature conduc- 
tance in terms of the effective S-matrix for low energy 
electrons via the Landauer formula. Below the Kondo 
temperature, a phase shift of 7r/2 occurs in the screening 
channel. To a good approximation, the low temperature 
S-matrix is simply determined by the unitary transfor- 
mation to the scattering basis and this n/2 phase shift 
in the screening channel. A small correction to this S- 
matrix occurs due to the Kr potential scattering term 
(or Vr in the case that Ed = —U/2). While this ap- 
proach confirms the results of Ref. [JJJ in the special case 
of a half-filled band in the leads, we find that changing 
the electron density in the leads has a large effect. 

We confirm some of these results by NRG calculations 
presented in § [Vj We only consider the simplest case in 
the Kondo regime, td = —U/2, symmetric coupling of the 
left and right leads to the dot, and with a half-filled band. 
We begin by completely describing the renormalization 
group flow of our model, predicting the form of the var- 
ious fixed points and crossover energy scales which are 
then confirmed in the NRG. Most notably, the Kondo 
temperature is extracted from the energy scale of the 
Wilson chain at which the crossover to the low energy 
strong-coupling fixed point occurs and agrees excellently 
with that predicted analytically. The effective S-matrix is 
compared to the low energy excitation spectrum over var- 
ious parameter ranges. We obtain quite good agreement 
through this comparison, including the small corrections 
from Vr. 



II. MODEL & ANALYSIS 

We start with a tight-binding model depicted in Fig.[T] 
The Hamiltonian for this model is 
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Each annihilation operator for the leads, cy, and the 
Anderson impurity, d, is a spinor where the spin in- 
dices are implied. The anti-commutation relationship is 

< c- , Cj> > — 5jj> . The number operator for dot electrons 

of spin /i is defined as n^ = didf,,. The various param- 
eters are described in Fig. [TJ We will assume that all of 
the couplings are real. The magnetic flux has been in- 
troduced through the parameter ip = 27r$/$o, $ being 
the magnetic flux threading the AB ring and $o = h/e 
being the magnetic flux quantum. We assume that the 
magnetic field generating the flux is small enough in the 



Sd 

□ 
t d -/ \t d+ 

t t t t .■' t t t t 

— O O O 6- - r -O O O O 

-4 -3 -2 -1 t 1 2 3 4 



FIG. 1. The lattice model described by the Hamiltonian of 
eq. GJ). 



vicinity of the wires so that we may neglect the Zeemen 
effect in the quantum dot and the leads. 

Although such a tight-binding model for the leads 
is not a very accurate description of leads in a semi- 
conductor heterostructure on which such geometries are 
often defined, we use it here as an example of a relatively 
simple model that contains a natural bandwidth of At. 

We now define a basis of even and odd combinations 
of electron operators 
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so that the Hamiltonian can be written as 
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where we have defined the shorthand notation 
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Hd remains unchanged. 

Immediately we notice that, for the case of zero flux, 
(p = 0, and symmetric coupling td- = td+, the model 
reduces to two decoupled chains, the even channel in- 
teracting with the quantum dot and having a potential 
scattering interaction —t' at j = 1 and the odd channel 
decoupled from the dot and with a potential scattering 
interaction t' at j = 1. However, in the general case of 
if 7^ 0, we must analyse both channels together. 

If we remove the dot from the system we are left with 
two independent channels, even and odd, with a poten- 
tial =Ft' at j = 1. As shown in Appendix [Bj this po- 
tential gives rise to two scattering phase shifts S k in the 
even/odd channel respectively, the form of which is given 
at the Fermi surface to be 
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where r' = t' ft and 5 
when kp = 7r/(2a), 5 + 



Note that, at half- filling 
= 6 where tan S = t' . 



These phase shifts will play an important part when we 
discuss the zero-temperature properties of this system in 

§nvi 

Noting that H |_ serves as a potential scattering term, 

we seek to transform to the scattering basis that essen- 
tially removes these interactions from the Hamiltonian. 
We do this by first introducing the complete set of wave- 
functions that solve the Schrodinger equation for Hq 
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and r\ is a positive, infinitesimal parameter. It is shown 
in detail in Appendix [X] that T kk , is given by 



■7± 
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(25) 

Thus, the Hamiltonian greatly simplifies in the q ek , q k 
basis to 



with a being the lattice spacing. We can expand our 
operators as 
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so that W k ,e k > \ = W k ,o k/ \ = S(k - k'). 
nian becomes 
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where we have defined 

Vkk> = ~ 



2t cos ka 
2t'a 



sinfca sin fc'a. 
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With H |_ written in such a simple form, we now trans- 
form to the scattering basis. Ignoring H t d for the mo- 
ment, we note that the only difference between the e and 
o channels is the sign of the v kk ' interaction. Hence, we 
define the scattering basis 
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The Hamilto- where 



#_+ = 
Htd = — 



dk- 



/ IV* r+*o t - 1* r^*J 

| ['tie 1 k U ek L do L k %k 



+h.C 



r 



dk' sin fc'a 



'„ ,A±(*0 



sinfca 



lir'e 



I a — iha 



(26) 

(27) 
d 

(28) 



(29) 



The last equality is proven in Appendix [X] 

With the potential scattering Hamiltonian H |_ van- 
ishing due to the transformation to the scattering basis, 
we are now free to rotate the basis once more to the chan- 
nel that couples directly to the impurity and its orthogo- 
nal complement. In this way, anticipating our discussion 
on the Kondo effect, we define the screening channel 
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where we have defined the asymmetry parameter 
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In this way the dot coupling Hamiltonian can be writ- 
ten as 
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This form of the hybridization was first found in ref. \l£ 



We are interested primarily in the Kondo effect which 



involves only the screening channel since it is the only 
one that couples to the quantum dot. Hence, one can 
perform the Schrieffer- Wolff transformatio n 2 ^ 23 so that 
the screening channel Hamiltonian assumes the form 



H= [' dke k V k CT ^f x 

Jo 
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dkdk' Jkk'Skk' ■ S d + K kk '^ 



kk'^k 



r V fc ?)(35) 



where we have defined the screening channel spin density 



Skk' 



and effective dot spin 
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with a being the three Pauli matrices (recall that \& 
and d are spinors). The coupling parameters are given 
by 



Jkk' = VdkVdk 

Kkk' 
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At this point, one can obtain a low-energy effective 
theory by integrating out high-energy modes in the usual 
way. The potential scattering term Kkk' is marginal and 
does not renormalize. We will discuss this term in more 
detail in § IIII Al and neglect it for now. The exchange 
interaction is relevant and diverges, giving rise to the 
usual Kondo screening of Sd by the screening channel 
Fermions for temperatures T below the Kondo tempera- 
ture Tk- There are, however, physical consequences due 
to the flux <p and inter-lead coupling t' that will be de- 
termined in § IIV1 

To summarize the analysis thus far, through a series 
of basis rotations we have cast the inter-lead Hamilto- 
nian into a potential scattering form. By transforming 
to the scattering basis, we have eliminated this potential 
scattering term and identified the operator that couples 
directly to the quantum dot. It is this combination that 

I 



will participate in the Kondo screening of the dot. Nev- 
ertheless, there are additional potential scattering terms 
that can arise in the screening channel and it is this sub- 
ject that we next discuss. 



III. ADDITIONAL POTENTIAL SCATTERING 

Our goal is to derive an effective theory of our system 
that is valid at low temperatures, keeping the leading 
order contributions in the effective strength of the Kondo 
coupling 



J 



1 



•// 
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(40) 
which has dimensions of energy and which we take to be 
a small parameter. 

The effective theory can be derived, to a first ap- 
proximation, by linearizing the dispersion Ek in a region 
— Q < k — kF < Q and approximating the coupling con- 
stants Jkk' and Kkk' by their values at the Fermi energy 
Jk F k F and Kk F k F - However, it will be shown that when 
the dot level is tuned to the value Ed = — U/2, Kk F k F 
vanishes to second order in V d k F ■ In this case, a more 
careful derivation of the low-energy Hamiltonian reveals 
that there is still an additional potential scattering gen- 
erated by the renormalization of Jkk' ■ This is higher 
order in J than the leading order contribution to Kk F k F 
written in eq. (I3TR) but contributes to the leading order 
term in the additional potential scattering when eq. (|39[) 
vanishes at Ed — —U/2, We address each of these cases 
separately below. 



A. Asymmetric dot Ed 7^ —V /2 

Restricting excitations to a small region about the 
Fermi energy as described above, the potential scattering 
term generated by the Schrieffer- Wolff transformation as- 
sumes the form 
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where 
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In order to observe the Kondo effect, we require that 
Ed ~ —U/2 so as to favor the formation of a local moment 
rather than a doubly occupied or unoccupied dot level. 
In this case, we see that Kr is of order J. However, 



for the precise value of Ed = —U/2, Kr vanishes and 
there is no potential scattering generated directly by the 
Schrieffer- Wolff transformation at low energies to linear 
order in J . 



The presence of this potential scattering term will give 
rise to a phase shift Sr at the Fermi surface in the screen- 
ing channel. As shown in Appendix [B] this is given by 



tan<5/? = —ttvKj 



(43) 



for small Kr and where v is the density of states at the 
Fermi energy. We will show in El UVBI how this additional 
potential scattering contributes to the T = conduc- 
tance of the AB ring. 



B. Symmetric dot Ed 



-U/2 



As discussed above, integrating out the high-energy 
modes to obtain a low-energy Hamiltonian leaves the 
marginal interaction Kkk' unchanged and so one obtains 
the term discussed in the above section. However, one 
can ask the question as to whether or not an additional 
potential scattering term is generated by the Kondo in- 
teraction Jkk' ■ Normally this is not the case for one often 
considers a Kondo interaction that is particle-hole sym- 
metric. It can be shown that this is not true for Jkk' of 
eq. (1381) . This is a consequence of a non-zero t' which 
necessarily breaks particle-hole symmetry. Although we 
have transformed away the explicit if interaction, the 
particle-hole symmetry breaking is manifest in this more 
complicated Kondo interaction. As a result, there is no 
symmetry forbidding this Kondo interaction from gener- 

I 



ating an additional potential scattering term and it is to 
the calculation of this that we now turn our attention. 

Consider a renormalization group scaling by integrat- 
ing out all of the wave vectors down to the Fermi energy. 
Although it is difficult to perform such a transformation 
exactly, one can make progress through a perturbativc 
expansion in J . The leading order contribution is of order 
J 2 which will be much smaller than Kr, eq. (|4"2"j) . which 
is of order J. However, Kr vanishes when Ed = —U/2 so 
that the J 2 term calculated below will contribute to the 
leading order term in the potential scattering. Hence, in 
this section, we assume that Ed = —U/2. 

Evaluating the Feynman diagrams to second order in 
the Kondo interaction Jkk' in eq. ([35]) , one finds a poten- 
tial scattering term generated of the form 



H R = V R 
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dkdk'% cr ^ k c , r 



(44) 



where the region of integration is restricted to small mo- 
mentum about the Fermi momentum kp and Vr is given 
by 



v '-h 



dk 
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The factor of 3/16 comes from the trace over spin de- 
grees of freedom and the denominator is simply the time- 
ordered propagator of the intermediate ^> s k CT Fermion. 



Substituting in the definition of Jkk' of eq. 
Vr can be written as 



38]) together with the definition of Vdk from eq. (|34|) and J from eq. ((40 
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The factors of I R are dimensionless integrals given by 
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with 
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To evaluate these integrals, we break them up into two regions 
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The imaginary parts from each integral cancel each other. Upon evaluation of the principle part of each integral, one 
obtains 
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Substituting this back into the above expression gives us our 

2>tt 2 {uJ) 2 1 + t' — 2jt' cos fc^a cos ip 
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final result for Vr: 
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where v is the density of states at the Fermi energy. Just 
as with the potential scattering term Kr, this Vr term 
will give rise to a phase shift in the screening channel 
given by 



tan Sr = —ttvVr 



(52) 



as shown in Appendix [Bl 

It should be noted that, although the potential scatter- 
ing generated by the Schrieffer- Wolff transformation van- 
ishes to order J, there may be a non-zero term at order 
J 2 in addition to that given by Vr calculated above. Such 
a calculation of the higher-order Schrieffer- Wolff terms is 
beyond the scope of this paper and so we leave it as a 
future project. 

In conclusion, the transformation analysis of section HI1 
provides a simple, generic way to account for the pres- 
ence of inter-lead coupling which takes the form of a po- 
tential scattering interaction. Such a transformation ef- 
fectively removes the potential scattering explicitly from 
the Hamiltonian in favor of a more complicated, particle- 
hole asymetric Kondo interaction when the dot is tuned 
to the Kondo regime. We have further shown that ad- 
ditional potential scattering terms are generated in the 
screening channel. The leading order contribution to this 
additional potential scattering is given by Kr, eq. (|42|) . 
in the case that Sd ^ — U/2 and by Vr, cq. (|5T]l . when 
£d = —U/2. In the next section, we analyse the physical 
consequences of this low-energy model. 



IV. PHYSICAL PROPERTIES 
A. Kondo Temperature 

One of the primary insights of the scattering transfor- 
mation analysis is in revealing how the Aharanov-Bohm 
ring influences the coupling between the quantum dot 
and the screening channel of electrons. That is, it allows 
us to obtain an expression for the dot-lead coupling in 
the Hamiltonian of eq. (|3"2"j) . given by Vdk F (hi the long 
wavelength limit), showing the dependence of the cou- 
pling on t' , ip, and kp. We then determine the t', if, and 
k F dependence of the effective Kondo coupling via the 
Schrieffer- Wolff transformation, eq. (1551) . This, in turn, 



r 



gives rise to a t', if, and kp dependent Kondo tempera- 
ture, the precise expression of which is easy to derive. 

Using the low-energy effective Hamiltonian, we deter- 
mine the effective Kondo coupling by evaluating eq 
at the Fermi energy 



J° = Jkwkp - V, 
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(53) 
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where J is defined in eq. (|40|) . The leading order RG 



definition of the Kondo temperatures* is 



T K = De- 1 ^"^ (55) 

Kondo temperature T$- = 



and dividing by the t' 
De -i/{iuJ)^ weget 



Tk t' 27 cos ip cos kpa + t' (l — 4 cos 2 kpa) + t 1 " 
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1 — 27T' cos ip cos kpa + t' 



(56) 

Although the denominator is always positive, we see that 
the Kondo temperature can be raised or lowered by the 
presence of the Aharanov-Bohm ring depending on the 
values of r' , (p and kp . This is shown in figures [5] and [3] 
which show the flux and t' dependence for various values 
of the other parameters. 

For the special case of half-filled leads, kp = ir/(2a), 
the result is particularly simple 



In — pt 



T 

2^J 



(57) 



showing that the Kondo temperature is independent of 
flux in this case. This limiting form of the Kondo tem- 
perature is verified by the NRG as discussed in § [VJ 



B. S-Matrix and Conductance 

The strong-coupling fixed point of the Aharanov-Bohm 
model under consideration can be described by a two- 
channel Fermi liquid. In this way, the fixed point model 
is fully described by a 2 x 2 S-matrix describing how the 
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FIG. 2. (Color online) The flux dependence of the Kondo 
temperature for a value of r' = 0.4 and 7 = 1. Here we see 
an increase in the flux dependence as the electron density is 
lowered. 



FIG. 3. (Color online) The r' dependence of the Kondo tem- 
perature for a value of Icf = 7r/6a. This exhibits the variety 
of behaviour that can be seen for different values of the flux 
and that the Kondo temperature always vanishes as t' -> 00. 



quasi-particle excitations of the two channels are scat- 
tered at the Fermi energy. In this section, we derive this 
S-matrix and relate it to the conductance between the 
two leads. 

The analysis of sections [TT| and Mil provide the fol- 
lowing simple picture of the strong-coupling fixed point. 
The direct coupling between the two leads, t' , gives rise 
to a phase shift (5 ± in the q e k and q k channels respec- 
tively. The form of these phase shifts is presented in 
eq. (|13p as computed in Appendix [B] By transform- 
ing to the scattering basis and removing the t' interac- 
tion from the Hamiltonian, we were able to identify the 
screening channel of eq. (|3TJ)) . Defining the orthogonal 
complement, ^| cr , to \&| cr and evaluating both at the 
Fermi energy (relevant here since we are talking about 
T = properties), we can write the relation between the 
screening/non-screening basis and the even odd basis in 
terms of the above phase shifts as 



Vj/scr 

o/scr 
fe 



= u 



Qek 

qok 



where 



(58) 






(59) 



with normalization 



N '= [(t 2 d _ +trf + )(sin 2 5 + (l+7Cos^)+sin 2 (5 (1 - 7 cos <£>))]' 



(60) 



I 

In the screening channel, there will be a phase shift phase shift. The second is the phase shift 5r generated by 
with two contributions. The first is the usual ir/2 Kondo the additional potential scattering, the leading order con- 
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tribution to which will either be Kr or Vrt^ depending 
on the value of Ed- Since the the additional potential scat- 
tering was obtained by integrating out the high-energy 
modes, the generated Hamiltonian term of eq. (|44|) must 
be considered as a low energy, long wavelength contin- 
uum model where the influence of the lattice is inconse- 
quential. The phase shift for such a model is derived in 
Appendix [B] and shown to be either that of eq. (l4"51 or 
eq. (Up. 

This is all of the information we require to write down 
the S-matrix in the even/odd basis: 



The far right matrix describes the potential scattering 
phase shifts due to t' in the q e k and q k channels, U 
rotates the basis to the screening channel and the matrix 
between U and U' describes the phase shift Sr due to Vr 
or Kr as well as the 7r/2 Kondo phase shift giving rise 
to the factor of — 1 = e 2l 5. 



Multiplying the matrices, we can write S as 



5' 



&oe &oo 



(62) 



S = U f 



„2l<5r, 



u 



J2io~ 





,-2ii5 



(61) 



with 



See = -Me 2iS+ (e 2lSn (l + ~fcos(p) sin 2 5+ - (1 - 7 cos ip) sin 2 <T ) 

Seo = -2Me i{s ~ +6++Sli \/3-ij sin ip) sin <T sin 5+ cos S R 

S e = -2Me i{5 ~ +s++SR) (j3 + j 7 sin <p) sin <T sin S + cos Sr 

S 00 = Me 2tS ~ ((I+7 cos p) sin 2 S + - e 2lSr < (1 - 7 cos tp) sin 2 6~) 



(63) 

(64) 
(65) 
(66) 



r 



where we have defined 



and 



t 2 



+ t 2 



(67) 



Next, we wish to form a combination of ip e and ip 
that corresponds to a right-moving wave incoming from 
the left. That is, we wish to form a superposition of 
the above two wave functions that has no left-moving 
component for x > 0. To this end, we form 



M 



(1 + 7 cost/?) sin S + + (1 — 7 cos v?) sin 6 



(68) 



To relate this S-matrix to the conductance, we first 
construct general scattering wave functions between the 
even and odd channels. Consider an incoming plane 
wave in the even channel that is then scattered into the 
even and odd outgoing channel according to the above 
S-matrix. Such a wave function takes the form 



^e = e- tklxl + Seee lk ^ + S oe sgn(x)e J 



(69) 



i> - o W'e ~ V'o) 



2 Wee + S 



(S e 



S e c 

-s n 



S 00 ) e %kx 

ikx 



,x > 
,x < 



(71) 



where, indeed, we find no e component in tp for x > 0. 

Looking at the x > portion of ^, we recognize the 
coefficient of the plane wave as the transmission proba- 
bility amplitude for transmission from the left lead to the 
right lead 



where the first term is the incoming wave in the even 
channel, the second term the scattered even wave and 
the last term the scattered odd wave. Similarly, consid- 
ering an incoming wave in the odd channel gives the wave 
function 

V'o = sgn(x)e- lk ^ + Se e lkM + S 00 sgn(x)e ik ^ . (70) 



J- — ~Z (See T S oe S eo S 00 ) . 



(72) 



Using the Landauer-Buttiker formula, we obtain an ex- 
pression for the conductance 



2e 2 
2e 



^f- { (I + 7 cos y) 2 sin 4 S + cos 2 ((5 
h 



5~ + Sr) + (1 - 7 cosy) 2 sin 4 (5~ cos 2 ((5 + - <5~ - <5r) 
sin 2 5~~ sin 2 (5 + [4 cos 2 Sr sin 2 93 - (l - 7 2 cos 2 ip) (cos [2((5 + - <5~)] + cos 25 R )] 
I [(1 + 7 cost/?) sin 2 <5 + + (1 — 7 cosy) sin 2 S~~\ 



(73) 



This is the most general expression for the conductance 
expressed in terms of the phase shifts 5 generated by the 
inter-lead coupling t' , the additional potential scattering 
Kr or Vr via Sr, and in terms of the flux tp. The latter in- 
cludes the explicit <p dependence written above as well as 
the dependence implicit in Sr via the flux dependence of 
Kr or Vr written in eqs. (l4"2l or (I5TJ) . Although the equa- 

I 



r 



tion is rather complicated, we see that the conductance 
satisfies the necessary symmetry relation G(ip) — G(—p). 
We now turn our attention to special limiting cases. 

For the case of kp — 7r/(2o) and td- = td+ considered 
in most previous studies, 8 + = —5~ = 5 with tan 8 = r' 
and the conductance simplifies to 



G\ 



2e^_ 
h 



cos 2 (2(5 - Sr) cos 4 ^ + cos 2 (2(5 + Sr) sin 4 ~ + cos 2 Sr snr 



2 



- (cos 4(5 
4 



cos 2(5^) sin 2 y 



(74) 
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FIG. 4. (Color online) The conductance as a function of dot 
level Ed for various values of the direct interlead transmission 
probability give by T5 = sin 2 25. Here, we assume the particle- 
hole symmetric value of half- filling, &f = it /2a and 7 = 1. 
This is the behaviour seen numerically in ref. IfOl . 



ft is interesting to compare this with the numerical re- 
sults of ref. [l(J. For the case of Ed 7^ —U/2, when Kr 
is the leading order contribution to Sr, we are able to 
qualitatively reproduce the Fano-Kondo behaviour seen 
in ref. [10| in the region Ed ~ —U/2 for which our analysis 
is valid. An example of this is given in figure [4] 

For the symmetric value Ed = —U/2 when Kr van- 
ishes, we can view the Sr generated by Vr as a small 
correction to the results of ref. LLCj. Indeed, in the limit 



(75) 



I 

of Sr — > and &f = 7r/(2a), our result reduces to 

GU^ = y(l-r 6 cosV) 

where Ti, — sin 2<5 is the transmission probability 
through the lower arm of the Aharanov-Bohm ring in 
the absence of the upper arm. This is precisely the form 
reported in [lOj for the case of a singly-occupied quantum 
dot. 

In this way, eq. (174)) can be viewed as an analytic de- 
scription of the results of ref. [f0|, the latter of which re- 
quired numerical input from the NRG. Such an analytic 
description is only valid for values of Ed close to —U/2 so 
as to strongly favor a local moment on the quantum dot 
whereas the results of ref. [l0| are valid for all Ed- On the 
other hand, our complete expression for the conductance, 
eq. (|73p . extends previous results to cases where the 
Fermi energy is not situated in a particle-hole symmetric 
manner relative to the band edges (e.g. kp 7^ it /(2a)) 
as well as taking into account the additional potential 
scattering Vr discussed in § IIIII 

To further examine the correction due to Vr, we look 
at the flux dependence of the conductance in figure [5] for 
the case that Ed — —U/2 and hence Vr contributes to the 
leading order behaviour of Sr. There, each of the differ- 
ent coloured lines indicates a different value of the direct 
inter-lead coupling t' as encoded by S. It is seen that the 
conductance contrast (the difference between the min- 
imum and maximum conductance) reaches a maximum 
for an intermediate value of the inter-lead coupling t' = 1 
(t'=t). 

Furthermore, it is shown that for r' < 1, the effect of 
the additional potential scattering Vr is to decrease the 
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FIG. 5. (Color online) The conductance plotted as a function 
of magnetic flux for &f = ir/2a, Ed = — U/2 and 7 = 1. Each 
of the different coloured lines indicates a different value of 
the inter-lead coupling r' . The solid lines are the prediction 
with Vr — (equivalently 8r = 0) with the dotted lines 
showing the finite Vr correction. This data assumes a value 
vJ = 0.287 for the bare Kondo coupling. 



FIG. 6. (Color online) The conductance plotted as a function 
of inter-lead coupling t p when kp = 7r/2a and Ed = —U/2. 
Here, J = 0.287 and <p = 0. 



conductance whereas for values r' > 1, the additional 
potential scattering serves to increase the conductance. 
This fact is made more evident in figure [5] where the 
conductance is plotted versus t' for <p = 0. There, one 
can clearly see the crossover from reduced to enhanced 
conductance around r' = 1. 

Given that Vr offers only a small correction, we look 
at the Sr = limit of the conductance for general kp 
which takes the form 



G|,5 H=0 = ^-[cos 2 (($ + -<r) (sin 4 <5+(l + 7 cos^) 2 + sin 4 <T(1 - 7 cos(^) 2 ) 
+ sin 2 <S~ sin 2 S + [4 sin 2 ip - 2cos 2 (<5 + - 5~) (l -7 2 cos 2 <p)] 
[sin 2 S + (l + jcosip) + sin 2 5~(1 — 7 cost/?)] 



(76) 



Even without including the small correction due to Vr, 
this is a generalization of the conductance reported in 
ref. [lfl which, like most similar studies, only considered 
the case where the leads exhibit particle-hole symmetry 
(kp = ir J '(2a) for our tight-binding leads). For quan- 
tum dots constructed on semi-conductor heterostructures 
where the two-dimensional electron gas has very low den- 
sity, the Fermi energy will be very close to the bottom 
of the energy band and so exhibit strong particle-hole 
asymmetry. Hence, the generalized forms for the con- 
ductance reported above seem to be more applicable to 
such devices than those reported in previous studies. 

The description of the conductance that emerges from 
this analysis is quite interesting. In the limit that r' — > 0, 
we recover the well-studied model of a single quantum dot 
embedded between two leads where one obtains unitary 
conductance at zero temperature. As one increases t', 



I 

interference effects play a stronger role until one obtains 
maximal interference at r' = 1 (TJ, = 1) where one is 
able to obtain total destructive interference in the form 
of zero conductance for certain values of the parameters 
(e.g. kp — ir /(2a) and ip = 0). As one further increases 
t', the transmission Tj, through the lower arm decreases 
and interference effects are diminished. 
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Fixed Point 



Free Orbital (FO) 
Local Moment (LM) 
Strong Coupling (SC) 



V d U Stability 



dot interaction 



Unstable 
oo Unstable 
oo oo Stable 



TABLE I. Summary of fixed points for the single-channel An- 
derson impurity model. 



V. SUPPORT FROM THE NUMERICAL 
RENORMALIZATION GROUP 

A. Phase Diagram &: Kondo Temperature 



1. Fixed points of the single- channel Anderson impurity 
model 



We begin by reviewing the various fixed points present 
in the single- channel Anderson model 26 before describ- 
ing the influence of the Aharanov-Bohm ring. When 
investigating low-energy, long-wavelength properties, it 
is customary to define a model in terms of continuous 
fields with a linearized dispersion relation characterized 
by a Fermi velocity vf- In this way, we can write the 
single-channel Anderson model in terms of right-moving 
ID electron annihilation operators ip{x) as 



H = v F I dxip^(x)(-id x )ip(x) + V d (ip\0)d + h.c.) 

(77) 



+ 2 (*<*-!) "2 



-C//2 (assumed 



where we have set the dot level to Ed 
throughout this section). 

This model has three fixed points summarized in Ta- 
ble HI The free orbital (FO) fixed point occurs when 
Vd = U = 0. This describes free ip Fermions with a de- 
coupled free dot level d. The spectrum of such a model 
is that of free Fermions plus the four degenerate, zero- 
energy states of the dot. 

The FO fixed point is unstable and flows towards the 
local moment (LM) fixed point as the energy scale is low- 
ered. The LM fixed point is characterized by a diverging 
U — > oo and Vd — 0. This LM fixed point is the same as 
the FO except that two of the four dot levels are ener- 
getically forbidden, namely, those for which cfd = and 
d/'d = 2. In other words, the quantum dot can only be 
singly occupied with either a spin up or spin down elec- 
tron. Hence, the spectrum will be that of free Fermions 
plus two degenerate, zero-energy states of the dot. 

The LM fixed point is also unstable and eventually 
flows to the strong coupling (SC) fixed point described 
by a diverging \Vd\ 2 /U — > oo. The nature of this fixed 
point can most easily be understood by first considering 
a Hamiltonian close to the LM fixed point with a small 
\Vd\ <C U. In this case, one can perform a Schrieffer- 
Wolff transformation^ 2 - perturbatively in Vd to obtain a 



H td + H dot ^J^{0)-^(0)-S d 



(78) 



where <x is a vector of the three Pauli matrices, Sd = 
d\cr/2)d is the effective spin of the singly-occupied dot 
level, and the coupling strength J is proportional to 
\Vd\ 2 /U. This is the Kondo interaction between the lo- 
calized spin of the quantum dot and the electrons in the 
leads. The SC fixed point of the Anderson model is es- 
sentially the same as the strong-coupling fixed point of 
the Kondo model wherein J — > oo and the dot spin is 
screened by forming a singlet with the lead electrons. 



2. Fixed points of the Aharanov-Bohm quantum dot model 

The low-energy transformations of section UU reveals 
that the renormalization group flow for the Aharanov- 
Bohm ring model under consideration will be very sim- 
ilar to that of the single-channel Anderson model just 
described. Indeed, we have learned that a single, inde- 
pendent combination of the lead electrons, ^ s ", couples 
directly to the dot just as in the single-channel Ander- 
son model. The precise nature of this screening channel 
will depend on both flux tp and the inter-lead coupling t' 
but the point is that there is a single channel available 
to screen the spin of the quantum dot. For simplicity, we 
consider only the symmetric case where t d - = td+ = t d 
and Ed = —U/2. 

The primary difference with the Anderson model dis- 
cussed in the previous section is the addition of some 
potential scattering phase shifts 6 depending on t' and 
the modification of the dot- lead coupling V d k F -^ Vdk F ■ 

We find that the FO and LM fixed points, with V d k F — 
0, will be the same as in the single-channel Anderson 
model with the addition of the phase shifts 5 arising 
from the direct tunneling between the leads which were 
incorporated into the definition of q e k and q k . The SC 
fixed point of the Aharanov-Bohm ring model will be one 
in which the dot spin is fully screened by the ^ s "' com- 
bination of lead electrons. Just as in the Kondo model, 
this will give rise to a ir/2 phase shift in the \I/| cr channel 
in addition to the phase shifts S ± arising from the direct 
tunneling t' . 

Furthermore, the FO and LM fixed points occur for 
Vdk F = and, since Vdk F encodes the t' dependence of 
the model, we predict that the cross-over scale of these 
fixed points will be unaffected by the presence of the 
Aharanov-Bohm ring {i.e. in the region of these fixed 
points, the t' and ip dependence of V d k F is inconsequen- 
tial). However, the cross-over energy scale to the SC fixed 
point, that is, the Kondo temperature Tk, will be influ- 
enced by the direct tunneling t' and flux p as discussed 
in§[IV3 

We can check these predictions for the fixed points 
of the Aharanov-Bohm ring model using the Numeri- 
cal Renormalization Group (NRG). This numerical al- 
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gorithm is exhaustively detailed in the pioneering pa- 
pers^— and in a recent review^ so we give only an 
outline sketch here. 

We begin with a long wavelength version of the Hamil- 
tonian described in eqs. (fl7]) - ([T9l) with a dispersion rela- 

I 



tion linearized about the Fermi energy ep = 0, a cutoff in 
momentum at k = ±Q (here, k is measured with respect 
to kp, i.e. k — kp — >• k), td- = td+ = td, and Ed = —17/2. 

The resulting Hamiltonian is 



H = vf I dkk \e\ek + o\oA -V p dk dk' [e\ek> - o\ok'j 

rQ 



+V d 



dk 



VJ 



Vj 



cos -e\ + sin -o\ ) d 



^(dU-iy 



(79) 



r 



where we have simplified our notation by defining 
the potential scattering term V p = —Vk F k F , Vd = 
—2td\/a/irsmkpa, and redehning the phase of o& so as 
to make all coefficients real. Note that this version of the 
Hamiltonian does not involve a transformation to scat- 
tering states. In this way, agreement between the NRG 
and results inferred from the transformations of sectionHIl 
will serve as support for the scattering transformation 
analysis. 

However, it should be observed that such a linear 
dispersion necessarily exhibits particle-hole symmetry 
whereas the tight-binding model discussed in section |TT] 
generally breaks particle-hole symmetry except for the 
special case of kp — it /(2a) that occurs when there is one 
electron per site. For this reason, the NRG as formulated 
here strictly serves only to support our scattering trans- 
formation analysis for the particle-hole symmetric case 
of kp — it /(2a). Nevertheless, we trust that our analytic 
results hold true for arbitrary kp. 

Setting up the NRG involves a series of transforma- 
tions and approximations that map the model for the lead 
electrons onto two semi-infinite tight-binding chains, of- 
ten termed Wilson chains, with hopping amplitudes that 
exponentially decrease with distance from the quantum 
dot 



TT 1 

n^(l + A" 1 )' E E A ~^« (/l/(n+i) b + h-c 

b—e,on=0 



D 




{Sd-lf-2uV v (flf Qe -flho 
(cos|/ t e + sm|/ t o )d + h. 



(80) 



In general, each Fermionic f ne and f no is a complicated 
linear combination of e^ and Ok respectively. The details 
of this relationship are not of great importance for the 
present discussion except to note that the Fermions cre- 
ated on the n — site by /J e and fl Q are proportional 
to the e and o electrons at the origin: /o e oc e(x = 0) 
and foo oc o(x = 0). Of the other parameters defined in 
this Wilson-chain Hamiltonian, 2D is the bandwidth and 
A > 1 is a dimensionless discretization parameter defined 



such that the continuum limit is recovered in the limit 
A — > 1. The dimensionless parameter £ n is given by 

e„ = (l-A-"- 1 )(l-A- 2 "- 1 )- 1 / 2 (l-A- 2 "- 3 )- 1 / 2 (81) 

and tends to unity for n >• 1. We have also defined 

T = 27wV d 2 (82) 

with v the density of states at the Fermi energy. 

The renormalization group is realized by truncating 
the infinite chain to N sites and rescaling the Hamilto- 
nian such that the eigenvalues are of order unity 



H N = A^- 1 '/ 2 



N-l 

EE 

b—e,o n—0 



A" 



kn (/„ 6 /(n+l)6 + kC. 



+U(d t d-1 



V P 



[jOeJOe JooJOo) 



C °S |/de+ sin l/oo)'' 



-h.C. 



(83) 



where the quantities with tildes are simply dimension- 
less versions of the original parameters of eq. (180[) with A 
dependent rescaling. The renormalization group trans- 
formation then takes the form of the recursion relation 

H N+1 = AJH N + £ N J2 (fLf(N+i)b + h.c.) (84) 

b— e,o 

and is realized by iterative diagonalization, using the 
eigenvalues and eigenvectors of Hn to define H^+i via 
eq. ([84)) . In practice, the eigenvalues are shifted so that 
the lowest one is zero. 

The finite Hamiltonian Hpf can be related to the 
Hamiltonian of eq. ([50)1 by 

Since the dimensionless scale of H^ is of order unity by 
definition, this indicates that the spectrum of Hn de- 
scribes the spectrum of the physical Hamiltonian at an 
energy scale given by 



E 



N 



Ul + A-^A-^-^^D. 



(86) 
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In this way, we can associate Hm with the effective 
Hamiltonian at the renormalization group energy scale 
En ■ Fixed points can be identified as regions of N over 
which the energy spectrum of the associated Hm changes 
very little (for unstable fixed points) or not at all (for 
stable fixed points). These fixed point NRG spectra can 
then be compared with that predicted by the scattering 
transformation analysis described above to test the va- 
lidity of said analysis. 



Our analysis of the fixed points follows that o f 26 ' 27 . 
Let us first consider the FO fixed point which, in terms 
of the NRG formalism, is defined by T — and U = 0, 

I 



resulting in 



AT-l 



s 1\ — J. 

H N , FO = A^- 1 )/ 2 £ J2 A-*£„ (flf {n+1)b + h.c 

^b— e,o n— 
~Vp yfoefoe - foofoo 



(87) 



This has the form of two decoupled Wilson chains, each 
with a potential scattering term at the origin. Such 
chains were analyzed ir*2I where the V p dependence of 
the single-particle energies was described in detail. 

Extending their analysis to two decoupled channels 
as described in eq. (IH7|) . one can diagonalize the non- 
interacting fixed point Hamiltonian and write it in terms 
of the single-particle and hole excitations 



H 



N,FO 



E &=e ,o Ei=t 1)/2 {vib(v P )gi b gnb + f n nb {%)h\ lh h n b 



J2 b 



Y.n=i y)i b (Vp)gi b 9nb + V nb (V P )hl b h nb j + f)+gl b g ob 



,7V odd 
, TV even. 



Here, g nb destroys a quasiparticle while h nb destroys a 
quasihole. The corresponding single particle/hole excita- 
tions are TV-dependent in general but, for N > 10 (ap- 
proximately) , they are found to only depend on whether 
N is even or odd, in which case one obtains ^-dependent 
energy levels V7 ± (V^,) or ?? ± (V^) respectively. 

The precise numerical values of these energy levels de- 
pend on A and V p . The V p dependence is described in 27 
where it was found that 



vUv P ) 



V~ b (-V P ) (89) 

and similarly for i). Furthermore, since the potential scat- 
tering in the e channel is equal in magnitude but opposite 
in sign to that in the o channel, the above relation can 
be written as 

vte(V P ) = vUVp) (90) 

and similarly for r). In this way, we recover a form of 
particle-hole symmetry even at finite V p where the energy 
spectrum of particles in the e channel are equivalent to 
the spectrum of holes in the o channel and vice versa. 

We can now combine these single-particle/hole excita- 
tions in multi-particle/hole combinations (being sure to 
respect the Pauli exclusion principle), together with the 
four degenerate zero-energy states of the dot level and 
so construct the FO fixed point spectrum. The lowest 
such energy levels are given in Table [H] along with the 
corresponding total charge Q and total spin S quantum 
numbers. 

The spectrum for the LM fixed point is closely related 
to that of the FO. The corresponding NRG Wilson-chain 
Hamiltonian for the LM fixed point is 



H 



7V,LM 



A (AT-l)/2 



N-l 



E 



n=0 



~Vp yfoehe - foofoo 

- lim U (dU - 1) 



(91) 
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which is identical to that for the FO fixed point with the 
addition of an infinite U Coulomb repulsion on the dot 
level. The corresponding spectrum of the LM fixed point 
will be the same as that for the FO fixed point with the 
exclusion of all of those states for which the dot level is 
empty or doubly-occupied as these now have an infinite 
energy cost. The lowest energy levels of the LM fixed 
point are listed in Tabic ITTT1 

To determine the spectrum of the SC fixed point, we 
must first identify the linear combination of electrons 
that screens the local moment on the quantum dot. How- 
ever, as discussed, we do not transform to scattering 
states in the NRG and so we simply use the combina- 
tion in eq. (|79[) that couples directly to the quantum dot 
as the screening channel, keeping the potential scatter- 
ing terms in the Hamiltonian, allowing the numerics to 
account for those terms directly. That is, we transform 
the original Hamiltonian, eq. (|79)). by rotating to a basis 



ip tp 

ipik = cos — e k + sin -o k 

ip2k = sin— e k -cos-o fe 



(92) 
(93) 



Yl A " 2 £" (flJ(n+i)b + h.c. J so that 
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Energy 


Num. Value Q 2S 


Energy 


Num. Value Q 2S 





0.0000 -1 


2??i + Via 


0.4485 1 




1 




1 




1 




1 2 

2 1 


vt 


0.1495 1 




1 


Vte + 2% 


1 




1 2 




-1 




2 1 




-1 2 


Via 


1 
-1 




-2 1 


2^+ + 27h~ 


0.5980 -1 




-1 2 




1 




-2 1 




1 


?vte 


0.2990 1 


*h + 


1.3580 1 




2 1 




1 




3 




1 2 


2r?r o 


-1 




2 1 




-2 1 


*?r e 


1 




-3 




-1 


Vle + Via 


-1 




-1 2 




-1 2 
1 




-2 1 








1 








3 








1 








1 2 







TABLE II. The lowest energies and associated total charge Q 
and total spin S quantum numbers of the free orbital (FO) 
NRG fixed point of the Aharanov-Bohm ring model for odd 
N. The numerical values for the single-particle excitation 
energies were obtained by diagonalizing the Hamiltonian of 
eq. (|87[1 using a value of V p — 3.0 and A = 2.5. All energies 
within a particular box are equal by eq. (|90[) . 



Energy 


Num. Value Q 2S 


Energy 


Num. Value Q 2S 





0.0000 1 


2^?i + % 


0.4485 1 
1 2 


nt 


0.1495 1 




1 2 


Vt + 2%~ 


-1 


■nio 


-1 
-1 2 




-1 2 


2i?£ + 2^ 


0.5980 1 


2^i 


0.2990 2 1 


??i 


1.3580 1 


2% 


-2 1 




1 2 


Vie + Via 


1 


Vie 


-1 




1 




-1 2 




3 







TABLE III. The lowest energies and associated total charge 
Q and total spin 5* quantum numbers of the local moment 
(LM) NRG fixed point of the Aharanov-Bohm ring model for 
iV odd. The single-particle energy levels are the same as in 
Table [TT] using V v = 3.0 and A = 2.5. All energies within a 
particular box are equal by eq. l]90p . 
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H = v F Y^ dkk ipl k ip bk + V d dk (ip\ k d + h 

'V fV'LV'iA' - i>\ k i>2k'\ + simp hpl k ip2k 



-Vp I dkdk' 



h.c. 



(94) 



and take V>ifc as the screening channel. 

The strong-coupling fixed point involves the ip\{x = 0) electrons forming a singlet with the dot local moment, 
effectively removing the V'i(O) and d degrees of freedom from the dynamics and giving rise to a n/2 phase shift in 
the ipi channel. One can then apply the standard NRG transformations and approximations to the resulting model 
in order to obtain a Wilson chain NRG form of the SC fixed point Hamiltonian. The ir/2 phase shift is implemented 
by shrinking the length of the ip\ Wilson chain by one site representing the removal of the site that is entangled in 
the Kondo singlet. 

The result is 



AT-2 



N-l 



ra=0 



H N , SC = M»-Wl J2 ^ {n+1)/ Hn (fUfn+1,1 +h.C.) + £ A_t en (fl, 2 fn + l,2 + h, 

-AsVp cos (pflxfai + Vp cos <pfl 2 fo,2 
-A-2 sjv p % sin (p (/ f ,i/o,2 + h.c. 



(95) 



Here, /„,i and f n ,2 are the NRG Wilson chain opera- 
tors derived from i\)\ and ipi respectively. The differing 
A prefactors are due to the normalizations required for 
the two different length chains. We have also added an 
additional factor, V' which arises from the additional po- 
tential scattering term in the screening channel discussed 
in section IIIII For now we simply take it as a single fit- 
ting parameter and return to its precise analysis in § IV Bl 
Since the ^2 channel does not participate in the screening 
of the quantum dot, we do not expect any additional po- 
tential scattering term proportional to /J 2/0,2- For the 
cross-term involving /g 1 /o,2 + h.c, we simply take the 
geometric mean of the two potential scattering terms of 



I 

the two channels and find that this provides a good fit to 
the NRG data. 



To obtain the SC fixed point spectrum, we first find the 
single-particle energy levels by numerically diagonalizing 
eq. (|9"5"|) for a finite value of N. As before, we find that for 
N > 10 (approximately), the energy levels depend only 
on the parity of N and not on its precise value. Unlike 
the FO and LM fixed point spectra, the resulting energy 
levels will depend on the flux ip in addition to the V p 
dependence. Similar to eq. (|88J). we can write the SC 
fixed point Hamiltonian in terms of the single particle 
and hole excitations 



H 



N,SC 



E 



(2iV+l)/2 
n=l 

(2JV+l)/2 






v t {V P , Vp, y)gig n + v n {V P , %, p)h\hA , A^ odd 

K{%, K, <P)gi9n + K{%, Vp, <p)hlh n ) , N even. 



(96) 



Because of the coupling of the 1 and 2 channels in 
eq. (|95|) . the quasiparticle excitations cannot be labelled 
by a channel index since it is no longer a good quantum 
number. 

The full many-body spectrum is constructed by com- 
bining these single-particle excitations in such a way as to 
respect Fermi statistics. The effect of the Kondo singlet, 
in addition to the 7r/2 phase shift already implemented 
in eq. (p?5"j). is simply to add an additional charge to the 
quantum numbers of the quasiparticle excitations due to 
the Fermion doing the screening. The lowest such ener- 
gies are listed in Table HVl 

Guided by the results of the transformations of sec- 



I 

tion [TT1 we have now identified the three fixed points 
of the Aharanov-Bohm quantum dot model and written 
the corresponding Hamiltonians in a Wilson chain form, 
eqs. (|5T|). (|9"Tj) . and (p?5"j). This allows us to determine the 
fixed point spectra, the lowest values of which have been 
listed in Tables |n] IIIII and IIVI We are now prepared to 
test these predictions by comparing these spectra with 
the actual energy levels that are computed in the NRG. 
This comparison is achieved by looking at the flow of 
the energy levels of each H^ (as defined in eq. (|83|l) for 
increasing N . An example is shown in Fig. [7] where we 
have plotted the lowest few energy levels of the Q = 1, 
5 = subspace as a function of odd N. It is shown 
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Num. Value 




Num. Value 


Energy 


H N ,sc NRG Q 2S 


Energy 


H N ,sc NRG Q 2S 





0.000 0.000 1 


l>2 + 2l/f 


0.8153 0.8158 -2 1 


v x 


0.0709 0.0711 1 


4 


0.8201 0.8203 2 1 


2v~ 


0.1416 0.1422 -1 


v+ + v x 


0.8910 0.8914 1 


v 2 


0.6737 0.6736 1 




0.8914 1 2 


V~ + l/f 


0.7445 0.7447 -1 


v+ + 2u~ 


0.9618 0.9625 1 




0.7447 -1 2 


2v- 


1.3474 1.3472 -1 



TABLE IV. The lowest energies and associated total charge Q 
and total spin S quantum numbers of the strong coupling (SC) 
NRG fixed point of the Aharanov-Bohm ring model for odd 
TV. The NRG parameters used are V v = 3.0 and <p — 1.047. 
The same parameters were used to determine the energy levels 
of Hn,sc where a value of Vp — 2.885 was found to reproduce 
the NRG data. 




FIG. 7. (Color online) The lowest energy levels with quantum 
numbers Q — 1, S = as produced by the NRG as a function 
of odd TV. The values of the predicted fixed point energies 
from Tables ITT1 and ITTT1 are indicated by arrows on the left side 
and energies from Table IIVI are indicated on the right. The 
parameters used to generate this plot are T/D = 0.0003142, 
U/D = 0.001, vV p = 1.05, and tp = 1.047. Here we see the 
unstable FO fixed point is approached for 5 < TV < 10, the 
unstable LM fixed point for 19 < TV < 33, and the stable SC 
fixed point for TV > 60. 



that the fixed point spectra predicted above are indeed 
approached in the appropriate regime. For example, for 
5 < TV < 10, all of the energies of the Q = 1, S = 
subspace of the unstable FO fixed point are approached 
with the proper numerical value as given in Table [TXJ 
Similarly, for 19 < TV < 33, the predicted energy levels of 
the LM fixed point (Tabic ITTT|) are approached. The same 
is true for the SC fixed point where, in Table HV1 Vp is fit 
in order to produce the fixed point spectrum produced by 
the NRG algorithm (for the parameters used to generate 
the NRG data, a value of V' — 2.885 was found to give 
the best fit). 

In Fig. [51 we show a similar plot of a single energy 
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FIG. 8. (Color online) The lowest non-zero energy level with 
quantum numbers Q = 1, S — as produced by the NRG as a 
function of odd TV. The different lines correspond to different 
values of the flux ip. The parameters used to generate this 
plot are the same as in Fig. [7) The value of TVk, related 
to the Kondo temperature via eq. (|100[) , is indicated by the 
arrow and is the same for all values of ip. 



level as a function of odd TV in the Q — 1, 5 = sub- 
space where the different lines indicate energies produced 
from different values of the flux p. Here wc sec that, as 
predicted, the FO and LM fixed point energy levels that 
are approached are independent of tp whereas those of the 
SC fixed point are strongly flux dependent. The slight 
flux dependence that appears in the LM region is proba- 
bly due to the fact that Vd is not quite zero {i.e. the LM 
fixed point is approached but never reached). Indeed, 
the flux dependence of the energy levels in this region 
decreases the closer the LM fixed point is approached. 

For a more quantitative analysis of this flux depen- 
dence, we plot the lowest NRG energy levels of the final, 
stable fixed point with those predicted by diagonalizing 
the Hamiltonian of eq. d95l as a function of ip in Fig. [9] 
The fact that a single parameter fit of Vp perfectly repro- 
duces the flux dependence of the entire NRG fixed point 
spectrum strongly supports the validity of the above RG 
analysis. Indeed, because the SC fixed point is stable, 
we can explicitly compare the fixed point spectrum pro- 
duced by the NRG with that predicted by eq. fl9"5*j) as we 
have done in Table ITVl for the first few levels. 



3. Kondo Temperature from the NRG 



In § IIV A[ we derived an expression for the Kondo tem- 
perature in terms of the inter-lead tunneling £', the flux tp, 
and the Fermi momentum kp, eq. (|56p . For the particle- 
hole symmetric value of kp — -k /2a, this expression takes 
the simple form of eq. ([57]) . It is this latter form that 
can be compared to the NRG which was derived from a 
model with a particle-hole symmetric linear dispersion. 

To do this, we must write eq. (|57p in terms of the An- 
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FIG. 9. (Color online) The lowest NRG energy levels of the 
final stable fixed point (red circles) are plotted as a function 
of flux ip and compared with those given by the fixed point 
Hamiltonian of eq. (|95[) (solid blue lines). The parameters 
used to generate this plot are the same as in Fig. [7J except 
here we use a value of vV v — 0.525. The single value of Vp 
was tuned in order to fit the NRG fixed point energy levels 
for ip = 0. This single parameter is able to reproduce the 
predicted flux dependence. 



derson model parameters appearing in the NRG Wilson 
chain form of the Hamiltonian, eq. ([5U|) . that serve as 
input to the NRG. First, we define an effective Kondo 
coupling for the continuum model of eq. (|79l) 



J 



2 M 
U 



(97) 



The factor of 2 is included because eq. ([79|) involves cou- 
pling to both the even and odd channels whereas, in 
eq. PU|). we defined J for the screening channel only. In 
transforming to the screening channel, a factor of v2 ap- 
pears in Vdk resulting in the derived J acquiring a factor 
of 2 which we account for here explicitly so that we can 
compare the NRG results with those derived analytically. 
Using the definition of T, (|8"2"j) . we get 



(98) 
2atr' /ir so that we 



vJ = 4T/(irU). 



Next, we recall that V p = —Vk F k f 

can write t' = ixvV v . The resulting expression is 



In — pt = 



2vJ 



(99) 



The right hand side of this equation, together with 
eq. (|98[) , now contains parameters related directly to the 
input parameters of the NRG. 

We now must extract the Kondo temperature from the 
NRG data for multiple values of V p and J in order to 
confirm the validity of (Ifjfll) . The Kondo temperature is 
defined as the energy scale at which the screening of the 
local moment takes place and the Hamiltonian crosses 



FIG. 10. (Color online) Select energy levels with quantum 
numbers Q — 1 , S — as produced by the NRG as a function 
of odd N. The different lines correspond to different values of 
the inter- lead coupling V p and the arrows of the same line type 
indicate the approximate value of Nk at which the Kondo 
crossover takes place for each case. The parameters used to 
generate this plot are the same as in Fig. [7] 



over to the stable SC fixed point. In the NRG, Tr- will 
be related to the value of TV at which the energy levels 
cross over from that of the LM or FO fixed point to those 
of the SC as described in the previous section. This value 
of N, which we denote Nk, at which the crossover takes 
place can be related to a corresponding energy scale using 
eq. (|5rj)l . namely 

k B T K « - (1 + A" 1 ) A-C^-WD. (100) 

One simply has to extract the value of Nk from the NRG 
energy level data in order to obtain Tk- In practice, we 
measure Nk for the lowest 20 NRG energy levels and use 
the mean value (Nk) to determine Tk- 

In Figure [TUJ we have plotted select NRG energy levels 
as a function of N for different values of V p . There is 
clearly a trend of increasing Nk with increasing V p which, 
from eq. (jlOOp . indicates a decrease in Tk as a function of 
V p as predicted in eq. (IM1) . Furthermore, if one looks at 
Figure [51 there is clearly no change in the value of Nk for 
the different values of flux ip indicating that there truly 
is no flux dependence in Tk when kp = it /(2a). 

For a more quantitive comparison, we have plotted the 
value of Tk extracted from the NRG as a function of vV p 
in Figure [Tl] for multiple values of the Kondo coupling 
J. The analytic form predicted in eq. (1991) provides an 
excellent, parameter-free fit to the numerical data. 



B. Phase Shifts and Vr 



As discussed in § IIVB1 the SC fixed point is comprised 
of two independent Fermi liquids characterized by two 
phase shifts. These phase shifts are determined by the 



18 






i-a-a-s_ai_ 


_ 


-0.5 




U 


1-. -1 












\o 




vJ = 0.400 
□ 0.255 


K 




0.127 
a 0.072 


No 


a -2 




< 0.051 
1 0.025 


\ 






Analytic Prediction 


\ "~ 


-2.5 






-3 


1 


1,1 



'1+ 



. K- 



FIG. 11. (Color online) A plot of Tk (symbols) as extracted 
from the NRG as a function of the input value of vV v . The 
different symbols describe data with different input param- 
eters giving rise to different effective Kondo couplings. The 
solid line indicates the prediction described in eq. (|99p . 



eigenvalues of the S-matrix of eq. (|62l) . In this section, 
we wish to compare these two predicted phase shifts with 
those derived from the NRG. 

Once again, given the particle-hole symmetric formu- 
lation of the NRG, we can only make this comparison at 
the special value of kp = 7r/(2a). In this special limit, 
one can see from eq. (TIB")) that 5 + = — 8~ = 8 where 
tan 5 — t' . We further simplify to the symmetric case 
td- = td+ = ta- In this case, the two eigenvalues of the 
5-matrix are 



+ channel 



channel 



with 



A± =~e i&R (iA±y/l-A 2 



A = cos 25 sin 6r + sin 2<5 cos 5 R cos ip. 



(101) 



(102) 



Writing these as pure phases A± = e 2la± , the phase shifts 
are given by 



cos 2a± = A sin S R qp V 1 — A 2 cos <5f 



(103) 



In the special case of ip = when the two channels fun- 
damentally decouple, one obtains 

cos 2a± = T cos(2(5 + S R ± S R ), {ip = 0) (104) 



or 



a_ 



\-6-6r, 


{<p = 0) 


(105) 


5, 


(p = 0). 


(106) 



The two phase shifts a± fully define the strong-coupling 
fixed point spectrum. 



FIG. 12. Energy level diagrams of the single-particle NRG 
energy levels of the two channels. The shift of each relative 
to the Fermi energy (here indicated by the dotted line) defines 
the phase shift in each channel. 



1. Phase shifts from the NRG 

First, we consider a system of two independent Fermi 
liquids on a finite line of length L and with linear dis- 
persion relations. The energy levels will then take the 
form 



r(0 



2nvp 
L 



(107) 



where q G Z and 6i are the phase shifts in the i th channel. 

The situation with the NRG is not quite so simple due 
to the non-uniform hopping in the Wilson chain that goes 
like hr n ' 2 at the n th site. However, one can still extract 
a sensible phase shift describing the overall shift of the 
(non-uniform) energy spectrum. We present a method for 
extracting these phase shifts from the NRG data that is 
similar to that used in 29 though we are much more mod- 
est about the claimed analogy between the non-uniform 
NRG spectrum and that of eq. (j!07[) . 

As discussed in § IV A 21 the many-body spectrum of 
the strong-coupling fixed point is built up of two channels 
of single-particle excitations, both of which we denoted 
together as z/„ where the ± superscripts indicate whether 
the excitation is that of a particle (+) or a hole (-). With 
knowledge of only the total charge Q and total spin S 
quantum numbers of each many-body energy from the 
NRG, one can identify the single-particle energy levels 
for each of the two channels which we denote as v nJr and 
v n _. It is from these that we estimate the phase shift in 
each channel. 

For clarity, let us assume that the NRG chain length 
N is even— and that the four lowest energy levels are 
ordered such that Vy, < v^_ < v^_ < i/y, , as depicted in 
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figure [12] The phase shift in each of the channels is going 
to be proportional to the lowest single-particle energy 
level in each channel, in this case, v^ + and t\_. However, 
because of the non-uniform A-dependent spacing of the 
energy levels, we normalize each phase shift by the lowest 
energy level spacing in their respective channels. That is, 
we define the phase shift as 



'i+ 



n_ 



(108) 



(109) 



If the channels are shifted in the other direction relative 
to the Fermi energy (that is, if the lowest single-particle 
excitation is that of a particle instead of a hole: 



'i+ 



< 



< I/, < !/,,), the phase shifts are taken to be 



a_ 



(110) 



(111) 



One can now extract the values of v u from the many- 
body NRG energy spectrum obtained by diagonalizing 
H jv as described in § IV A 21 Assume that N is sufficiently 
high such that the RG has reached the strong-coupling 
fixed point. The ground state, describing no particles or 
holes and set arbitrarily to Eq =0, will have total spin 
quantum number S — and a charge quantum num- 
ber of Q = +1 or Qo — — 1 depending on whether the 
lowest single-particle energy is a hole or a particle respec- 
tively. Let us assume that Qo = +1 for clarity. Then, the 
values of v^, and v^_ are given by the two lowest many- 
body energies with a charge quantum number of Q = 
and spin quantum number S = 1/2 (these lowest ener- 
gies would be Vy + and z/jt_ if Qo = — !)• The values of 
r/J, and v^_ are given by the lowest many-body energies 
with charge quantum number Q — +2 and spin quantum 
number 5 = 1/2 (in the case of Qo = — T ^f+ ari( ^ v \- 
would be given by the lowest Q = —2, S = 1/2 many- 
body energies). In this way, one can extract the single- 
particle/hole energies and estimate the phase shifts from 
the NRG data. 

As an illustration of the A dependence of these phase 
shifts, we consider the simple case of zero-flux, (p = 0. In 
this case, the original Hamiltonian can be completely de- 
coupled into two separate channels and so the two chan- 
nels operate completely independently. The channel cou- 
pled to the quantum dot is the screening channel and 
so obtains a 7r/2 phase shift in addition to that given 
by 5 whereas the other channel is non-interacting with 
only a potential scattering phase shift —6. This can be 
seen clearly in figure [T3] where we have plotted the two 
phase shifts as a function of V v (recall from § IV A 21 that 
V p — —Vk F k F and so is related to t' via eq. (J2"TJ0 . 

The most striking feature of figure Q~3] is the different 
A dependence in the phase shift of the screening and 




FIG. 13. (Color online) The phase shift of the two channels 



as a function of V v 



-Vkpk, 



from the NRG data with an ef- 



fective Kondo coupling v J = 0.255 and zero flux. In this case, 
the two channels are independent with the screening channel 
phase shift starting at 7r/2 for vV v = and the non-screening 
channel phase shift starting at for vV p = 0. The symbols 
are the phase shifts obtained from the full NRG many-body 
energy levels. The ascending solid lines are the phase shifts 
obtained from the single-particle energy levels of a single non- 
interacting Wilson chain with potential scattering V p and the 
descending solid lines are 7r/2 minus the ascending lines. The 
solid lines do not take into account the small correction due 
to the additional potential scattering Vr that occurs in the 
screening channel. 



non-screening channels obtained from the NRG data. To 
help understand this, we have plotted as solid lines the 
phase shifts that one would expect in a non-screening and 
screening Wilson chain (we ignore the effects of the small 
correction due to Vr for now). For the non-screening 
channel, one can diagonalize directly the Wilson chain 
Hamiltonian with a potential scattering V p at the first 
site using different values of A and so obtain the single- 
particle energy spectra directly without having to per- 
form the NRG. From this direct single-particle spectra 
one can define the phase shift as described above and 
these are plotted as the solid ascending lines. As can be 
seen, these match perfectly the phase shifts in the non- 
screening channel obtained from the NRG data, as they 
must. 

To leading order (again, neglecting Vr), one might ex- 
pect the phase shift in the screening channel to be sim- 
ply 7r/2 minus the above A-dependent phase shifts since 
the potential scattering in the screening channel is equal 
in magnitude but opposite in sign to that in the non- 
screening channel. We have plotted this expectation as 
the descending solid lines in the figure. On the contrary, 
the phase shifts obtained directly from the NRG data 
show very little A dependence compared with the non- 
screening channel. The precise reason for this is unknown 
though it may be due to a similar A dependence in the 
additional potential scattering Vr as shown in figure [TS] 
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FIG. 14. (Color online) The phase shifts of the two channels 
as a function of V p . The symbols denote phase shifts derived 
from the NRG data while the lines are the analytic predic- 
tions. The solid black line is the curve expected if there is 
no additional potential scattering (i.e. Vr = 0). We have set 
ifi = to generate this plot. 



that is compensating for the A dependence of the bare 
5 phase shift. Despite this, we nevertheless obtain good 
support for our prediction of the phase shifts from the 
tight-binding model. 

In the remainder of our analysis, we will use the A 
dependent phase shift obtained from diagonalizing the 
potential scattering Wilson chain discussed above for 
the bare phase shift 5 generated by V p that appears in 
cqs. (fTU2"|) , (|T05|) . and (fTUo) . SeeSI for more information 
on the V p dependence of the NRG spectrum. 



2. NRG evidence for Vr 

We now turn our attention back to the additional po- 
tential scattering Vr that was derived in section lHTl Hav- 
ing shown that the phase shifts can be extracted from the 
NRG, we can now compare the predicted phase shifts in 
eq. (|103|) with those of the NRG. For simplicity, we con- 
tinue to assume t c i- — td+ — id and Ed = —U/2. To 
compare our analytic results with those of the NRG, we 
use the same correspondence as was used in section lV A 31 
namely r' = irvV p and vJ — 4T/(wU). 

We focus first on the case of zero flux, ip = 0, where 
the phase shifts take an especially simple form given in 
eqs. (|105|) and (I106p . These two phase shifts are plotted 
in figure [14] as a function of V p where the symbols indi- 
cate those values derived from the NRG data while the 
lines are the analytic prediction from the tight-binding 
model. Here we see that, indeed, only the phase shift 
of the screening channel (the one that obtains tt/2 when 
V p = 0) deviates from the Vr = prediction, indicating 
that an additional phase shift is generated in the screen- 
ing channel only. However, Vr provides only a small cor- 



FIG. 15. (Color online) The additional potential scattering 
Vr as derived from the NRG phase shifts (symbols) and from 
the analytic tight-binding model (lines) for various values of 
the effective Kondo coupling J. We have set ip — to generate 
this plot. 



rection so it is easier to extract Vr from the NRG phase 
shifts and compare its functional form directly with that 
of eq. floTi 

To extract Vr, we take the arctan of the derived NRG 
phase shift and subtract from that the tt/2 contribution 
arising from the Kondo screening as well as the bare 
phase shift S due to V p . This latter phase shift will be 
A dependent and can be calculated numerically as de- 
scribed in21. 

In figure [151 we compare directly the predicted depen- 
dence of Vr on Vp with that derived from the NRG phase 
shifts for various values of J. We find that both analytic 
and numeric calculations of Vr share the same qualitative 
behaviour, peaking around vV p sw 0.3 (corresponding to 
t' ~ t in the original tight-binding model), but that pre- 
cise quantitative agreement is not obtained. The nature 
of this disagreement is discussed in section [VTl 

We next turn our attention to testing the (vJ) 2 de- 
pendence in eq. (|5Tj) by plotting the value of Vr as deter- 
mined from the NRG phase shifts versus v J in figure 1161 
The most striking characteristic is the apparent deviation 
from universal behaviour as U/D approaches unity. We 
see that this trend is captured by the U dependence in 
eq. (|51l) but that precise quantitative agreement is elu- 
sive, perhaps because of the presence of a cubic term 
which we do not consider. A complete analysis of Vr 
with an Anderson impurity rather than reducing, via the 
Schrieffer- Wolff transformation, to one with a spin impu- 
rity may elucidate the nature of this behaviour. 

For further analysis, we fit the largest data set with 
U/D = 0.001 to a third degree polynomial of the form 



Vr = a + ai(vJ) + a 2 (vJ) 2 + a 3 (ujy 



(112) 



A third degree polynomial was chosen instead of a sec- 
ond degree function because the data goes to quite large 
values of vJ where we expect our second order analysis 
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FIG. 16. (Color online) The additional potential scattering 
as determined from the NRG phase shifts as a function of 
effective Kondo coupling J. Each different symbol uses a 
fixed value of the dot Coulomb repulsion U while varying V 
such that the range of J, given in eq. (|98p. is roughly the same 
for each iteration. The solid line presents the best fit third 
degree polynomial to the U/D — 0.001 points. In this data, 
0. 



Coefficient 


Value 


ao 


-0.00081 


Ol 


0.045 


«2 


-0.62 


«3 


0.50 


&2 


-0.42 


h 


0.29 



TABLE V. The parameters for the best fit of eqs. (fTT2|) 
and (fTT3)) to the U/D = 0.001 data in figure [T6l 



FIG. 17. (Color online) The phase shifts of the strong- 
coupling fixed point as determined from the NRG (symbols) 
and compared with that predicted in eq. (|103[l (lines). Here, 
the effective Kondo coupling is v. J = 0.191 and vV v = 0.25. 
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to break down. The values of the parameters are tabu- 
lated in Table [V] It is seen that the coefficients ao and 
Oi, which we predict to vanish, are indeed at least an 
order of magnitude lower than the quadratic and cubic 
coefficients. Doing another fit neglecting these first two 
terms, that is, to a form 



V R = b 2 {iyJ) 2 +b 3 (isjf 



(113) 



gives b-2 = —0.42 which is the same order of magnitude 
as the value of —0.24 predicted by eq. (J5l"j) . 

Up until this point we have been focussing primarily 
on the form of the additional potential scattering Vr and 
so, for simplicity, have taken the flux ip = 0. In figure [T71 
we have plotted the phase shifts a± versus the flux (p as 
derived from the NRG with comparison to the predicted 
form described in eq. (|103l) . There we find the agreement 
to be quite good and suggests that our predicted flux 
dependence is robust. 

Finally, we note that, although it seems that the A 
dependence of the screening channel phase shift is sup- 
pressed (see figure [T5]) , there does appear to be some 



FIG. 18. (Color online) The value of the additional potential 
scattering Vr as derived from the NRG for various values of 
A, all using a value of J — 0.510. 



systematic A dependence in Vr itself as seen in figure IT51 
In order to take this effect into account, one would need 
to derive an expression for Vr from the Wilson chain 
Hamiltonian, eq. (|80l) . as opposed to the much simpler 
tight-binding model as was done in eq. (|51[) . We still find 
convincing agreement of the behaviour of Vr between 
that derived analytically and from the NRG despite this 
apparent A dependence. 



VI. DISCUSSION 

We have presented a systematic study of a minimal 
model of an Aharanov-Bohm ring with an embedded 
quantum dot connected to two conducting leads. Al- 
though aspects of such a system have been studied by 
other groups in the paslj^r— , our work provides a com- 
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plete picture of the physics of such a system when the 
quantum dot chemical potential q is near — U/2 and the 
system is in the Kondo regime including new effects not 
discussed previously. 

In particular, we have elucidated precisely how the 
Kondo effect arises in such a system by identifying the 
screening channel; we have completely mapped out the 
renormalization group flow of the system and its depen- 
dence on flux ip and inter-lead tunneling £'; we have cal- 
culated the dependence of the Kondo temperature and 
conductance on the same parameters as well as, for the 
first time, the electron density in the leads (via the fac- 
tors of kp appearing throughout); we have calculated the 
effects of additional potential scattering that arises from 
the breaking of particle- hole symmetry; we have provided 
wide numerical support from the NRG for many of our 
findings that goes beyond simply computing the occu- 
pancy of the quantum dot as ir>i£ or the dot density of 
states^. Although our work is quantitatively precise, 
the physical picture that arises has been stated in simple 
physical terms that fully describes the zero-temperature 
properties. 

It is interesting to compare our results for the Kondo 
temperature with those of refs. [13| and [lj. In the for- 
mer reference, the authors use a slave boson mean field 
theory to estimate the Kondo temperature for variable 
sized rings. For the smallest configuration with only one 
site in the ring in addition to the quantum dot, they find 
a flux dependent Kondo temperature assuming particle- 
hole symmetric leads, kp — it /(2a). Although the cal- 
culation of ref. [l3| was for a different model than that 
considered here, the two models are quite close and the 
nature of this apparent discrepancy is not clear. It is 
interesting to note that the authors of ref. [l3| find very 
similar behaviour at kp = 7r/(2a) to that found by us 
for electron densities less than half-filled, kp < ir/(2a) 
(see figure^]). It may be that the particle-hole symmetry 
breaking caused by moving away from half-filling in our 
calculation mimics the particle-hole symmetry breaking 
caused by the negative on-site energy of the additional 
site in the ring used in ref. [lj- Perhaps it is this type 
of particle-hole symmetry breaking that leads to a flux 
dependent Kondo temperature. This is speculation and 
further analysis of both methods would be required to 
resolve this apparent discrepency. 

Reference [lj follows a very similar procedure to that 
used here, transforming to the scattering basis and iden- 
tifying the screening channel. However, they mainly 
consider the U —¥ oo limit with finite dot energy level 
Ed- Their subsequent scaling analysis, assuming half- 
filled leads with particle-hole symmetric Fermi energy 
Sp = 0, produces a flux dependent Kondo temperature. 
Although this seems to contradict our conclusion that 
the Kondo temperature is flux independent at half-filling, 
our result was obtained in a much different limit, with 
Ed ~ —U/2. The authors do claim that, for finite U, the 
flux dependence is suppressed (though still present) when 
Ed = —U/2. However, we find no evidence of any flux de- 



pendence in the Kondo temperature when kp = it /(2a). 

We close our discussion with a few comments on the 
apparent discrepancies presented in the NRG evidence 
for the additional potential scattering Vr. As discussed 
in the text, we expect there to be cubic and higher order 
contributions to Vr that we do not calculate so discrepen- 
cies for values of v J that approach unity should be ex- 
pected. However, discrepancies remain even for relatively 
small values of v J and we offer here some possibilities for 
why this might be. 

As written at the end of section IVB21 the correspon- 
dence between the tight-binding model used to derive Vr 
in eq. (|5ip and that used in the NRG is only approximate, 
especially for values of A > I. This leads to artificial A 
dependence in many of the quantities extracted from the 
NRG as has been presented above. This is probably true 
for the value of Vr extracted from the NRG, as seen in 
figure HS1 suggesting that the form of Vr may be non- 
universal in that it may depend on the details of the 
band structure of the leads. 

To explore the universality of the form of Vr, the au- 
thors have repeated the derivation of Vr for a model with 
a linear dispersion in the leads rather than the tight- 
binding cosine dispersion presented in the text. It was 
found that, while qualitatively the same as the form of 
Vr in eq. (|5ip . the two forms of Vr did differ in numeri- 
cal details. From this we conclude that the form of Vr is 
non-universal. In light of this fact, one would ideally re- 
peat the calculation of Vr, not for the tight-binding chain 
presented but for the full A-dependent Wilson chain and 
so obtain the A dependence of Vr. However, given the 
non-uniform 'tunneling amplitudes' in the Wilson chain 
that go as A~ rl ' 2 for hopping from the n th site, such a 
calculation would be very difficult. 

Another possible source for this discrepancy is the pos- 
sibility of additional contributions to potential scatter- 
ing arising from the Schrieffer- Wolff transformation. We 
have performed such a transformation to second order in 
Vdk and concluded that the potential scattering Kk F k F 
that arises vanishes when Ed = —U/2 so that, in this 
regime, Vr contributes to the leading order term in the 
potential scattering. However, given the fact that a non- 
zero t' breaks particle-hole symmetry, there is nothing 
preventing the Schrieffer- Wolff transformation from gen- 
erating a potential scattering term that is fourth order 
in Vdk (equivalently, second order in J). It would be 
interesting though non-trivial to carry out the Schrieffer- 
Wolff transformation to higher orders to see if indeed 
such potential scattering terms are present and if they 
can account for the disagreement with the NRG. 

Despite all of these possibilities, it is clear that such 
a Vr term is present in both the tight-binding model 
as well as in the NRG and that they share the same 
qualitative behaviour and modestly agree quantitatively. 
Given this, we expect such a Vr term to be present in 
any real physical system and we expect it to share the 
same qualitative dependence on flux (p, inter-lead tun- 
neling t' (peaking around t' w t) and on electron density 
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in the leads via the dependence on Uf but do not claim 
that it will be precisely as that given in eq. (UTTl) which is 
based on an overly simplified tight-binding model. Fur- 
thermore, although present, the contribution of Vr to 
the conductance is very small for typical values of v J, as 
seen in figure [5j and so will probably be difficult to de- 
tect explicitly in any physical system. Nevertheless, the 
remainder of our analysis is robust and confirmed numer- 
ically and provides a framework in which to think about 
such quantum dot systems. 
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Appendix A: Details of the transformation to the 
scattering basis 

Let us consider the Hamiltonian H = Hq + H |_ where 

Hq and H h are given by eq. (fTT|) and (fT5)l respectively. 

We will now demonstrate that, under the transformation 
of eqs. (|2"2"|) - (|2l?)) . the above Hamiltonian takes the form 



H = I dk e k (q\ k q ak + q\ k qbk 



(Al) 



Another way of saying this is that we demand the trans- 
formation to be such that 



H,qL 



£ kq' 



,a = a,b. 



Substituting into eq. (|A2[) the definition of q* , in terms 



of el and o[ from eqs. (l2"2"j) - (j2U]) and using the relations 



H,e\ 



H,ot 



e k e\ 



dk'v k ' k e\, 



ii 



e k o\ - I dk'v k 'kO k , 
Jo 



(A3) 
(A4) 

(A5) 



one obtains 

Jo 

Substituting now the definition of (/>,, from eq. ([24]) in 
the above expression gives an integral equation for T kk , 

qk 




If we now take the ansatz 

Tr,, = Trj sin ka 



(A7) 



and substitute this into eq. (JA6I) together with the def- 
inition of v kk i from eq. (1211) and e k — — 2t cos ka, one 
obtains the following equation for T k 



m+ It' a sin ka 

T h = T ; — 



(A8) 



where, as before, r' = t' /t and I k is the dimensionless 
integral 



4 = - 



dy 



sin y 



cos y — cos ka + 177 



(A9) 



This integral can be solved in the complex plane. Making 
the change of variables z = e ly , one can write this as 



h 



1 



dz 



(z^-iy 



z 2 (z 2 — 1z (cos ka — irj) + 1) 



(A10) 



where the contour of integration is the unit circle centered 
at the origin in the complex z plane. 

The integrand has poles at z = zq = (second order) 
and at z — z± (simple), the latter given by 

; 1 -i • 1 1 ■ cos ka . . 

z± = cos ka ± 1 suite — ir\ =F V~i r- (All) 

|sinKa| 

Since r\ is a positive infinitesimal quantity, one can show 
that \z+\ < 1 whereas \z~\ > 1 so that only the z + and 
zq poles lie within the contour. Applying the residue 
theorem 

h = -~ (Res(z = z Q ) + Res(z = z+)) (A12) 



7re 



-ika 



(A13) 

Substituting this final value back into eq. (|A8|) and the 
resulting T k back into the ansatz, eq. (|A7I) . produces the 



(A2) promised form of T t 



kk' 

m+ 2t'a sin ka sin k'a 

T? u = =F- 



- kk' 



■k 1 ± r'e 



/ —ik' ■ 



(A14) 



as stated in eq. (|23|). 

Finally, we compute the functional form of T k as de- 
fined and stated in eq. ([2l?]) . Substituting the definition 

as defined in eq. (|2"4"|) into the definition of Tt 



of 



A±W 



and using the derived form of T kk ,, eq. (|A14[) . gives 



1 k — 



dk' sin k'a </> fc , 



= sin ka I 1 =F 



7r 1 ± r'e" 



-ika 



(A15) 
(A16) 



where I k is the same dimensionless integral defined in 
eq. (IA9|) and computed in eq. (|A13|) . Hence, the final 
result is obtained upon substitution 

sinfca (A17) 



rj = 



lir'e 



t o—ika 



as reported in eq. 
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Appendix B: Potential scattering phase shift 



Continuum Model 



1. Lattice Model 

Consider a single semi-infinite tight-binding chain with 
an on-site potential at the first site 



H=-t 



j=i 



0* 



■j+i 



- t'e[ ei 



(Bl) 



This is the same as the even-channel Hamiltonian H = 
H + H-+ of eqs. fl?| ) -(fl8 ]l in the limit t d+ = td- = 0. 
The presence of a finite t' will give rise to a phase shift in 
the single-particle wave function and it is the calculation 
of this phase shift that is the subject of this appendix. 
We write the eigenvectors of the Hamiltonian as 



At low energies (long wavelengths), one can take the 
continuum limit of the tight-binding model and linearize 
the dispersion relation about kp. In this way, one can 
write an approximate real-space Hamiltonian 



H = v F dx^{x){-id x )^{x) + V R V t (0)V'(0) (B9) 



where vp is the Fermi velocity. We assume that vVr <C 1 
where v is the density of states at the Fermi energy. 

As in the lattice model, we introduce eigenvectors of 
the Hamiltonian 



\<h 



dxMx)^(x)\0) 



(BIO) 



.7=1 



10} 



(B2) which satisfy the Schrodinger equation 



where {<fij} are coefficients to be determined such that 
they satisfy the Schrodinger equation 



HI 



£k\ 



(B3) 



with £fc = — 2t cos ka. The Schrodinger equation can be 
written as the following series of algebraic equations 



t(j) 2 - t'tfii 



£k<pl 



-t( 



D,-_l 



+ <Aj+i) = e k<f>j 



,J>1. 



(B4) 
(B5) 



The equations in the second line can be solved by taking 
coefficients of the form 



4>j = sin (kja + 5k) 
and 8k is determined by eq. (|B4|) to be 

t' sin ka 



tan<5i, = 



(B6) 



(B7) 



1 — r' cos ka 
where, as before, r' =t'/t. This is the form of the phase 



shift St that occurs in the even channel. The phase shift 



in the odd channel is the same but with r' — > — r' so that 

t' sin ka 



t&nSi 



±- 



T =F r' cos ka 



(B8) 



H\^ k )=v F k\k). (Bll) 

This puts the following condition on the functions </>& (x) 

- iv F d x (j)k{x) + V R (j)k(0)5(x) = v F k(pk{x). (B12) 



We now take the ansatz 





| 


e i(kx+S R ) 


,x>0 


(t>k(%) 


" 


e i(kx-5 B ) 


,x < 




1 


cosSr 


,x = 


ith derivative 








| 


ike 


,i(kx+S R ) 


,x > 


d x (f>k(x) = < 


ike 


,i(kx — 6n) 


,x < 


1 


2ismSjid(x) 


,x = 0. 



(B13) 



(B14) 



In order for <j>k(x) to satisfy eq. (|B12j) . we require 

Vr 



tan<5r 



2v F 



= —-KvVn 



(B15) 
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